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ABSTRACT 

The evolution of multiple stellar systems can be driven by Kozai cycles and tidal friction (KCTF), which shrink the orbit of the 
inner binary. There is an interesting possibility that two close binaries on a common long-period orbit experience mutually-induced 
KCTF. We present the discovery of a possible new quadruple system composed of two unresolved eclipsing binaries (EBs), CzeV343 
(V~ 13.5 mag). We obtained photometric observations of CzeV343 that completely cover the two orbital periods and we successfully 
model the light curves as the sum of two detached EBs. We provide confidence intervals for the model parameters and minima timings 
by bootstrap resampling of our data. One of the EBs shows a distinctly eccentric orbit with a total eccentricity of about 0.18. The two 
orbital periods, 1.20937 and 0.80693 days, are within 0.1% of a 3:2 ratio. We speculate that this might be the result of KCTF-driven 
evolution of a quadruple system and we discuss this hypothesis in the context of other quadruple systems composed of two EBs. We 
make our double EB fitting code publicly available to provide a tool for long-term monitoring of the mutual orbit in such systems. 

Key words, binaries: close - binaries: eclipsing 



1. Introduction 

The stability of multiple stellar sy stems requires that the stars 
are hierarchically organ ized (e.g. lEggleton & Kiseleval 1 1 9951: 
ISterzik & Durisenl U998) with an inner binary orbited by one 
or more outer bodies, although our knowledge of the forma- 
tion and evolution of these systems is still limited. In the case 
of hierarchical triple systems, the presence of an outer body can 
drive the evolution of the inner binary throug h the Kozai cycles 
and tidal friction (KCTF) mechanism (e.g. iMa zeh & Shaham 
| 1979HKiseleva et alJl998hlEggleton & Kiseleva-Eggletonl200U 
iFabrvckv & Tremainel2007l) . Kozai cycles periodically raise the 
eccentricity of the inner binary while the tidal friction efficiently 
dissipates the orbital energy during the close pericenter pas- 
sages. The net result is a decrease of the orbital period of the 
inner binary. The KCTF can in principle explain the existence of 
close binaries with orbits smaller than their pre-main-sequence 
dimensions. Indeed, for orbital periods P < 5 days the frac- 
tion of binary systems with a distant companion appears to 
increase (e.g. Tokovinin etlril 120061: |P ribulla & Rucinski 2006; 
ID' AngeloetalJ 120061: iRucinski et alJ 12007) and the eccentric- 
ity of the inner orbit is generally low (e.g. | Duquennov & Mavor 
ll99lHRaghavan et al.ll20iol:lDong et all2012ah . 

One prediction of the KCTF mechanism is the existence of 
a distant companion on typically high-inclination orbit around 
the inner binary. The parameters of the outer-body orbit can 
be obtained eithe r directly via astrom etry or spectroscopy (e.g. 
iHorn et al.1 119961: ISoderhielml 1 19991) . by using the inner bi- 
naryas a clock that orbits around a common barycenter ( e.g. 
lMaverlll990l: iBorkovits & Hegeduslll99fl iGies et aill2012l) or 
by co mbination of both (e.g. Rib as et alj|2002h [Zasche & Wolj 
2007). Furthermore, the properties of the inner binary (such 
as the apparent inclination) can change due to the perturba- 
tions from the outer bod y, potentially leading even to a cessa- 
tion of the eclipses ("e.g. lEggleton & Kiseleva-Eggl etonll200lt 



IZasche&Paschkd 120121) . In all cases, an often heterogeneous 
set of data has to be analyzed, which requires correct estimation 
and appropriate treatment of observational uncertainties. 

One of the possible arrangements of multiple systems are 
quadruple systems composed of two close binaries on a mu- 
tual long-period orbit. There is an interesting possibility that 
the two close binaries mutually influence their orbits through 
KCTF, although to our knowledg e no theoretical studi es of 
such systems have been performed. iTokovinin et al.l (12003) sug- 
gested KCTF as an origin of high eccentricity in 41 Dra, which 
forms a quadruple system with 40 Dra. A potential wealth of 
information on these systems is available if both binaries ex- 
hibit eclipses. To date, only three such system s have been con- 
firme d as SB4 binar i es: BV Dra and BW Dra dBatten & Hardiel 
[19651 iBatten & Lul fl986l). V 994 Her dLee et alJ 120081) and 
KIC 4247791 dLehmann et alll2012l) . The first two are visual 
bina r ies wh ile the third is currently unresolved. Additionally, 
lOfirl (12008b proposed that OGLE J051343. 14-691837.1 is also 
an SB4 binary where the two binary periods are in an ex- 
act 3:2 ratio, but later investigations revealed only an SB2 
with unexplained spectral changes related to the se cond pe- 
riod dKolaczkowski et al.l l20ld iRivinius etaf] 1201 ll) . OGLE- 
LMC-ECL- 16549 is c omposed of two unr esolved binaries with 
very different periods dGraczv k et al. 201 1 ). There are two other 
quadruple systems with only one of the binaries favorably ori- 
ented to show eclipses (Sh kolnik et al.l 120081: iHarmanec et al.l 
120071) . 

In this paper, we present the discovery of CzeV343, a double 
eclipsing binary similar to V994 Her and KIC 4247791, but with 
the periods of the two eclipsing binaries very close to a 3:2 ratio. 
We describe CzeV343 and our photometric data in Section [2] In 
Section [3] we present a model of the observed light curves and 
in Section|4]we discuss the physical connection between the two 
binaries. In Section|5] we discuss and summarize our results. 
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2. Observations 

The photometric variability of a V ~ 13.5 mag star CzeV343[3 
was detected during a search for new variable stars with 0.25 m 
//5.4 Newtonian telescope equipped with G4- 16000 CCD cam- 
era^ and a coma-corrector located at a private observatory in the 
Czech Republic. Initially, no photometric filter was used to max- 
imize the throughput because the camera is sensitive between 
350 and lOOOnm with a maximum at 550 nm, but later we ob- 
tained V and Ic band photometry with a higher quantum effi- 
ciency G2-3200 CCD camera. Most of the exposures were 180 s. 
All images were calibrated with appropriate dark frames and flat 
fields created as a median of five individual dark and flat ex- 
posures. We performed differential aperture photometry using 
C-Munipaclfl which is based on DAOPHOT dStetson|[T987l) . A 
nearby star of similar brightness and similar color (GSC 02405- 
01305) was chosen as a comparison star to minimize the effects 
of differential extinction. In total, we obtained 833 useful photo- 
metric measurements in 14 nights spanning ~ 80 days before the 
star became unobservable as it moved into conjuction with the 
Sun. Our data are summarized in Table [1] 

CzeV343 is a previously unknown variable star. Initially, 
CzeV343 appeared as a typical detached eclipsing binary in our 
data. However, subsequent observations revealed that the light 
curve is highly peculiar with three types of minima, as can be 
seen in the phased ~ 1 .2 day period light curve shown in Figurefl] 
During the course of our observations, the minima at phases 
0.25, 0.6 and 0.9 were slowly drifting with respect to the pri- 
mary and secondary minima. The only explanation that fits the 
observed periods and the light curve is that CzeV343 is com- 
posed of two eclipsing binaries with orbital periods of 1 .2 and 
0.8 days, similar to V994 Her and KIC 4247791. 

3. Model of the light curves 

In this section, we discuss the modeling of our data. The ob- 
served data consist of magnitudes m, and their uncertainties cr, 
measured at times f, distributed in M datasets. Each dataset con- 
sists of Nj measurements from one night and in one passband. 
We converted f, to the barycentric Julian dates in barycentric dy- 
namic al time (BJDt db) using the on-line t oofl of Eastman et al.l 
d2010). Similarly as Lehmann et ail (|2012|) . we modeled the ob- 
servations as a sum of the fluxes of two detached eclipsing bina- 
ries, 

u 

m(ti) = -2.51og w [F A (a A ,ti) + J3F B (a B ,ti) + Fe] + Yj c Afi>W 

,/=i 

where Fa and Fb are the fluxes from eclipsing binaries A and 
B scaled by the factor /3, and Fq represents any additional flux 
in the system. Each eclipsing binary is described by a set of pa- 
rameters a. For each binary we varied the orbital period P, time 
of primary minimum To, the sum of radii r\ + measured rel- 
ative to the semi-major axis, the ratio of radii r2/Vi, the surface 
brightness ratio cr, the inclination i of each orbit with respect 
to the observer, and for the binary with longer period we also 
vary two eccentricity parameters (e sin co, e cos a>). We fixed the 
mass ratios of both systems to unity, because this has little effect 
for detached systems. Lacking the spectral types of the compo- 
nents and given the quality of the photometry, we set the gravity 



1 GSC 02405-01886; ff=5 h 48 ra 24!008, cS=+30°57'03'.'64 (J2000) 

2 Parameters of our detectors are given at http://www.gxccd.com. 

3 http://c-munipack.sourceforge.net/ 

4 http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html 
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Fig. 1. Light curve of CzeV343 folded with the orbital period 
of system A of about 1.209 days. Each dataset is plotted with a 
different symbol explained in the right part of the plot. 
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Fig. 2. Light curves of the two eclipsing binaries with the model 
contribution of the other component removed phased with the 
appropriate P and 7q. The best-fit model (Table|2| is shown with 
solid red lines. The dashed red line shows best-fitting model with 
eA = 0. The smaller panels show the magnitude residuals phased 
with respect to the period of each of the systems. 



and linear limb darkening coefficients for all components and all 
filters to 0.5. We allowed for magnitude shifts cj between the in- 
dividual datasets and the model; Aji is unity if the time instant f, 
is within the dataset j, and zero otherwise. 

We fit all parameters simu ltaneously by minimiz i ng y 2 
using the routine cmp£it0 (iMorel Il978t iMarkwardtl [2009) 
with FA( aA,tj) and FB(an,tj) computed with the code 
JKTEBOP dPopper & Etzellll98U ISoufhworth et aTll2004l 120071 
iBruntt et aLll2006l) . We performed a number of bootstrap resam- 
plings of the original data to derive reliable confidence intervals 



http://cow.physics.wisc.edu/~craigm/idl/cmpfit.html 
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Fig. 3. Comparison of SDSS DR7 images of CzeV343 (left 
panel) and a similarly bright star about 1' away (right panel) 
with isophotal contours overplotted in blue. We show images in 
the z' band, where CzeV343 does not saturate and the brightness 
difference between the two binary systems should be the small- 
est. 

of the parameters. For both the original and bootstrapped data 
we started the minimization at a number of random positions in 
the parameter space to find the lowest^ 2 . 

We found that our model has degeneracies in several param- 
eters. First, the scaling factor /3 is to a large extent degenerate 
with the additional flux Fc and several of the binary parame- 
ters such as inclination. We performed the minimization for a 
range of fixed values of f3 with either Fc free to vary or fixed 
at Fc = 0. With Fc left to vary, we found quite small Ay 2 as 5 
for 0.6 < ft < 1.8. Because there is little prior information on /3, 
we decided to fix Fc = 0, which yields f3 & 1 .85 with a second 
minimum at ft ~ 1.50 differing by A^ 2 * 10. The assumption 
of equal fluxes ft = 1 is Ay 2 « 40 worse than ft ~ 1.85. In our 
fiducial solution, we set Fc = and ft = 1.85. There is also a 
degeneracy in r2/r\ in the sense that bootstrapping results occa- 
sionally yield two peaks: one with r2/n < 1 and the other with 
rilr\ > 1. The relative probability of the two peaks depends 
on ft. This degeneracy occurrs, because r^jrx essentially sets the 
flatness of the primary and secondary minima, which are con- 
strained only by a few datapoints, and on the amount of blended 
light of the other binary, which is controlled by ft. For/? = 1.85, 
about 86% of the bootstrap resamplings have > 1, which 
we then used to obtain confidence intervals on the parameters. 
Our best-fit model has x 1 = 1040 for 833 measurements in 16 
datasets and 30 free parameters. The model parameters are given 
in Table |2] and the light curves of the two eclipsing binaries are 
shown in Figure|2] The model is a good fit of the data. 

4. Physical connection between the binaries 

Even though a combination of two eclipsing binaries fits the 
data well, this does not prove that there is a physical connec- 
tion - the two binaries can simply be projected on the same 
position on the sky. To see whether CzeV343 can be resolved 
with our data, we analyzed centroid positions relative to a 
nearby star as a function of the flux. This method is used to 
const rain false positives in the search for tr ansiting exoplanets 
(e.g. lBatalha et alJ2010tlJenkins et all2010h . However, the plate 
scale of our setup is l'.'39/pixel with a typical stellar FWHM of 
3 pixels. Furthermore, our setup was designed for a wide field 
of view rather than for precise astrometry. As a result, we can 
only limit any shifts to be smaller t han about ff.'S. Additio nally, 
we analyzed the SDSS DR7 images dAbazaiian et alj |2009, plate 
scale of 0.396"/pixel and a typical stellar FWHM of ~ 2 pixels) 



that have a resolution superior to our data. In Figure [3] we show 
CzeV343 in comparison with a close star of similar brightness. 
There is no discernible difference between the two stars or any 
other stars nearby. More specifically, the measured PSF elliptic - 
ities and position angles are the same. CzeV343 is thus unre- 
solved given the available data. 

5. Discussions and conclusions 

We showed that the light curve of CzeV343 is well described by 
a sum of the fluxes of two eclipsing binaries, but it is unclear 
whether the two eclipsin g binaries are physi cally related. In the 
case of KIC 424779 l. lLehmann et alJ ( 120 12h argued for a phys- 
ical connection based on the similarity of masses, spectral types 
and periods of the four stars. Similar reasoning can be applied to 
CzeV343 as well, because the temperatures of the four stars must 
be similar given the lack of systematic differences between the 
V and Ic band residuals (Fig. 0. However, there is an observa- 
tional bias against discovering blended double eclipsing binaries 
with significantly different temperatures without any regard to 
the physical connection between the components. 

The question of whether systems like CzeV343 and 
KIC 4247791 are truly quadruple systems can be established 
by obtaining accurate distances to all four stars, which requires 
precise photometry, radial velocities and spectral classifications. 
Alternatively, the physical connection can be proved by observ- 
ing the effects of the gravitational interaction between the two 
systems with the added advantage of determining the parame- 
ters of the mutual orbit. In particular, light travel time effects in 
both systems with the same period would determine the mutual 
orbit. This method requires precise minimum timings with real- 
istic uncertainties and with proper separation of the flux of both 
components. For example, there are over 50 minima timings of 
V994 Her over the last ~ 20 yearfl but their practical useful- 
ness is very low because in some cases simultaneously occurring 
minima were not properly separated. To start the long-term mon- 
itoring of CzeV343, we determined O-C values for our datasets 
by fixing the binary parameters to the values in Table [2] and fit- 
ting for magnitude shifts cj and changes in To- We bootstrap- 
resampled data in each dataset and repeated the fitting to give 
reliable 90% confidence intervals. The results are presented in 
Table Q] Because fitting such unusual light curves is not com- 
monly done, we make our code publicly available^ with the hope 
that this will allow determination of minima timings of Cze V343 
and similar stars in the future. 

The physical parameters of the two binaries (Table |2) are 
not extraordinary in any way and are to some extent degenerate 
given our data, but it is worth pointing out that system A has an 
eccentric orbit. The secondary minima occur at phase ~ 0.52, 
but they also last longer than the primary minima, which yields 
a total eccentricity of e& « 0.18. The best solution with = 
is worse by Ay 2 w 1360 and is shown in Figure|2]with a dashed 
red line. An eclipsing binary with e « 0.18 and P ~ 1.2 days is 
quite rare (iDevorl 12005). Observations of the apsidal motion in 
system A and orbital pre cession can constrain the parameters of 
the mutual orbit as well (Eg gleton & Kisele va-Eggleton 2001). 

Finally, the most interesting property of CzeV343 is that the 
two orbital periods are very close to a 3:2 ratio: 

IP 

1 - — - = (8.5 ± 0.2) x 10~ 4 . (2) 
j "b 



6 http://var.astro.cz/ocgate 

7 http://www.astronomy.ohio-state.edu/~pejcha/czev343 
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Table 1. Summary of observations. 
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Notes. Each line corresponds to a single dataset. <fl) Beginning of ob- 
servations in BJDxdb _ 2 450 000. {b) Duration of observations in hours. 
(c) Passband of the dataset. C stands for unfiltered observations while V 
and l c are standard Johnson-Cousins filters. The fitted minus mean 
value of To in the units of 10~ 3 days. The uncertainties are 90% confi- 
dence intervals. w Primary minimum. Secondary minimum. 



The difference is visible in Figure [TJ as a slow drift of system B 
minima. Given their relatively small orbital separations, systems 
such as CzeV343 likely experienced KCTF that shrunk their or- 
bits. There is a possibility that during this process the orbital pe- 
riods became locked in a 3:2 resonance. This interpretation is ap- 
pealing because the light curve of OGLE J051343. 14-691837.1 
also appears as a combination of two contact binaries with pe- 
riods in an exact 3:2 ratio, alt hough this has not been con- 
firmed yet through spectro scopy (IOfirl2 008 ; Kolaczk owski et all 
2010: [Rivinius et alJ[20TTh . The periods in V994 Her, 2.08 and 
1.42 days, are considerably farther from 3:2 with a calculation 
similar to Equation © giving ~ 2 x 10~ 2 . For KIC 4247791, 
the periods of 4.10 and 4.05 days are clearly off the 3:2 ratio, 
although their relative difference is again ~ 10 -2 . The nearly 
resonant periods in CzeV343 may very well be a coincidence, 
especially if there is no physical connection between the bina- 
ries. There is no theoretical basis for such a behavior neither, 
but to our knowledge, there has not been any study of KCTF- 
driven evolution of a quadruple system. How ever, with the cur- 
rent interest in astrophysical Kozai cycles (e.g. Blaes et al. 2002; 
Fabrvckv & Tremaine 2007; Karas & Subr 2007; Perets & Naoz 
2009t lThompsonll201 li iDong et alj|2012bl) expanding the stud- 
ies to binary pairs may be a logical extension. 
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Table 2. Parameters of CzeV343. 
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2450000. {b} Fixed during fitting. 
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